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ABSTRACT 



The results of coherence measurements of three free precession 
proton magnetometers are discussed. The measuring and data process 
sing instrumentation is described. Ihe noise associated vri.th the 
instrumentation has been measured and presented in statistical fom. 
The results of measurements with elementary local perturbations are 
included to illustrate the target detection problm. The limitations 
of the instrument for target detection are discussed and imporvements 
suggested. 

The writer wishes to express appreciation to Professor Carl E. 
Menneken for suggesting the problems, and Professor Harold A. Titus 
for encouragement and assistance in the preparation of this thesis. 
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1, Introduction 



The use of magnetometers for target detection requires identic 
fication of signals in the presence of the noise associated with the 
earth* s magnetic field. The signal, in mangr cases, will be smaller 
in magnitude than the normal diurnal variation of the field, local 
fixed perturbations and geomagnetic piilsations. Ihe system noise 
is of major concern in a detection scheme using magnetometers. This 
paper presents the preliminary results of a stucfy of three proton 
free precession magnetometers operated under non laboratory condi=> 
tions at Naval Postgraduate School La Mesa Village site. 

The free precession magnetometer consists of a coil of wire 
in which a source of protons is placed. A strong field is applied 
to orientate the proton spins in a direction approximately nomal 
to the earth's field. The polarizing field is then removed in a man= 
ner such that the spins cannot follow the field. The earth's field 
then causes the protons to process about the lines of magnetic force 
at the Lamour frequency given by 2 




where '^p s gyro magnetic ratio 

s magnetic raoment/angular momentum (1) 



and 

|-j = magnitude of external field 

(earth's field) 

The net precessing magnetic moment causes a voltage to be induced in 
a coil surrounding the sample. Since the processional frequency is 
dependent only on the external magnetic field, the measurement of H 
consists of measuring the frequency of the induced emf. The inherent 
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accuracy of the free precession magnetometer is determined by the 
spectral line width (l)s which is related to the relaxation time 
T 2 > and is given bys 



where 



L 







= 3.t 

H s Half intensity width 
T 2 s 1 sec (meas) 
s 10°^ GAUSS 



( 2 ) 



The gyromagnetic ratio has been raeasxired by the NBS to better than 
1 part in 10^, Full advantage of this inherent accuracy has not been 
utilized because of the problems associated with the frequency measure- 
ment. The frequency measurement is complicated by the fact that only 
a limited time is available for measurement. The signal received 
from the magnetometer is in the few microvolt rangej hence locally 
generated electronic noise and transients from the polarizing opera- 
tion cause the signal to noise ratio to be small. Standard counting 
techniques are used for measurement. To get the required accuracyj 
the frequency must be obtained from period or multiple period average 
measurements. Ihe process is digital in natxire and lends itself 
very nicely to computer processing of the data. 

The main disadvantage is that any reasonable attempt to obtain 
an analog signal proportional to the measurement is an analog of 
period rather than frequency? therefore, the analog is the inverse 
of the quantity of interest. This does not present serious problems 
for target detection since the absolute magnitude of field is not 
of interest. 
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2. Ins trumenta tion 



2,1 Measurement 

The engineering aspect of the measurement problem can be stated 
simply as follows? 

Given ? Three transient type signals of frequency approximately two 
kc/s occuring at the same time 5 with usable duration of about 3/U 
sec. The signals are in the microvolt range. The recurrance rate 
is to be 2 sec. The measuirements must be conducted in the presense 
of high electronic noise such as relay transients^ contact noise 
and transients from the collapsing polarizing field. 

Find ? The frequency of the oscillations accurate to 1 part in 10^ 
or better, and suitable means of data aquisition and processing. 

The funds available excluded the possibility of using commer^ 
cially constructed instruments or components. Also, modifications 
and interfacing wovild have to be done if components were used, Hsnce, 
the basic measuring instrumentation was constructed at NPGS, No 
detailed circuit design is presented| however the circuit diagrams 
are included in Appendix IV, 

Assume the precession signal has been amplified in a suitable, 
low noise, narrow band preamplifier, and converted to a signal hav= 
ing fast waveform. Referring to the simplified block diagram of 
Fig, 1, the basic counter operation is as follows? Let the slow 
and fast counters consist of flip=flop counting units, connected to 
form a binary counter. The gate control flip=flop is capable of 
being triggered into either state. The multiplier converts the clock 
frequency to a harmonic of the precision 100 kc/s oscillator. 
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BASIC COUMTER UNIT 



Input 




100 kc 

Ctysta.1 Oscillator 



Figttre 1 

Let all counters be reset and ready for operation. The unknown 
input signal is applied to the gate. After a delay of approximately 
1^0 msj an external signal triggers the gate control to open the sig=* 
nal gate. The long delay is necessary to allow for damping the polari= 
zing field and relay transients. The first pulse through the signal 
gate opens the clock gate. The fast counter then counts the number 
of clock pulses occuring during the time interval determined by the 
slow counter. The slow counter is increased one for each cycle of 
xinknown input signal. After the slow counter has reached a predeter= 
mined number, both signal and clock gates are closed, and the niamber 
in the fast counter is a value of the tine interval of 2® cycles of 
unknown input signal, or can be reduced to the average period of the 
signal. 



k 



Let us define 

N s number in fast counter 

2*^ s number of cycles of unknown signal 

tp g period of clock (^^sec.) 

tjj g period of unknown signal 



Then we have 



N irf. 




= N tp jz"' 

Substituting in (1) gives 

H = 

Yp l^/^p 

and with n s 10, 



(3) 




2345 4oo X 

N 



Ox) 



The accuracy of such a scheme depends ong 

a. The accuracy in determining the zero crossing of the unknown 
signal. 



b. Propagation delay in the slow counter. 

c. Gating speed. 

d« Accuracy of the clock. 

e. Plus or minus the las bit, due to noncoherent gating as 
illustrated in Fig. 2. 
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EFFECT OF NONCOHERENT GATING 
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Figure_2 



The errors due to determining the zero crossing are usually the domi^ 
nant error in using this techniquei however when averaged over 102k 
cycles they become small and the largest error is the plus or minus 
last bit, as shown. Ihe accuracy in determining the zero crossing 
in the presence of noise (2) isg 



where? 

t^ g maximum error in trigger time in sec. 

Vn s peak to peak noise in volts. 

S s slope of ideal signal at point of trigger levels 
For a sinusoidal inputs 



V. 



n 




s 





Assume a signal to noise ratio 2 Vg/7n s 20 db, and 



w ^ !c c /s 



then? 
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For 102ii periods 8 

O. 3 X /O ^seo 

The plus or minus last bit error if the clock is 100 kCg is ^10 
sec., and the other errors mentioned are negligible compared to the 
10^,,4^sec« error. 

Consider the bit limitation in detail for the earth’s field 
50?50 ^ , wtiich is the approxinate value at Konterey, California, 
lets 

N s OOI 3 U 05 I 3 ^ 1 bit 
1 bit 10 ^^,^sec. 

them 

/i H gii.i y 

This becomes the sensitivity of the insiruraentc If the clock is 
increased to 1 me /sec, 1 bit s ands 

^ H e 0.11 'y' I however other noise becomes appreciable and 
cannot be neglected. The circuits were designed for 1 m/see opera 
tion, but the measurements to date vrere conducted at 100 kc/s and 
300 kc/s clock rates. The added complication is to construct 3 
meters that id-ll read identical, the input circtiitry' and slow count 
ers must have identical characteristics. This is Impcssible 
a practical point of view. The I'es’alt is that Uie three counters 
Kith an ideal signal input can vary-^ 1 bit in comp-arirjg their read 
ings. 

A simplified block diagram and photo of the nagnetometer is 
shown in Figs, 3, 4 and 5. It must be remerfbered t’nal operation 
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is to be in remote locations where extremes of temperatures line 
voltage fluctuationss poor gro\mdss etc* will be found| hence cir= 
cuits that are sensitive to these variations must be carefully com® 
pensated* 

For the experiments under considerations locating the sensing 
coils approximately 500* apart would be sufficientj however for fur® 
ther experiments a distance of greater than 500* might be necessary. 
Therefore, a provision to place the preamps in the line nearer the 
sensing coils for improved signal to noise ratio was included. Bias 
for the preamps is applied to the signal line, from the console to 
eliminate the need for battery replacement. The relays that provide 
the two to three amps of polarizing current, and provide for dissapa® 
tion of the stored energy when polarizing current is turned offs must 
precede the preamp. Since this unit is placed near the sense coil, 
it must be magnetically clean. 

Referring to the block diagram the components that are common 
to the three units are the clock, timer, recorders and the power 
supplies. 
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COUNTER INSTRUlffiNT 
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2.2 Data Aquisitton and Processir^ 

rw m ~p ~r^i I »I ~ « nr m r\irm rriw^iMT^ i^limiii fr 

Figvire 6 shows a block diagram of the data collection and repro- 
duction system. 

The scheme used for data handling is certainly not optimum in 
tlie sense of tape utilization and ease of operation^ but an Ampex 
CPIOO \ms available and required no additional expendit-ure for tape 
recorder. Also vdth minor modifications a data, recorder can be uti- 
lized at a later date if funds become available. 

Considei' the form of the data. The sensing coils are polarized 
for one second, at the anl of this time the polarizing current is 
turned off, the transients are damped| and the signal is fed to tlie 
countii^ units. After 102U cycles of the input signal (about 1/2 
sec.) the measurement is complete. In the fast counters, there are 
three binary nxmibers to be processed. The magnitude of the binary 
numbers, for H s 50950 vd.th clock frequency s 100 kc/s, is about 
0013lilU6g or l6 bits. With the clock incx'eased to 300 kc/s, the 
number is 00U25260g or 18 bits. Hence, at the end of the counting 
interval, there are 1x8 bits to be recorded in a format foi' computer- 
input. The Avtpex CPIOO is a "portable" instrumentation recorder 
xd.th lU parallel inputs and capable of speeds from 1 7/8 ips to 60 
ips. All of the record and reproduce amplifiers that were on hand 
wei*e not capable of recor'ding D.C, levelsj therefore, the information 
was put on tape in the form of pulses. Also, for reasonable spsod 
in reproduction, the data was recorded at 1 ?/6 ips and repr-oduced 
at 60 ips. A record of approximately 6,5 hours can be reproduced in 
10 minutes. The machine used to read the raw data is a CDCI60 com- 
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DATA AQUISITION AND PROCES3IN3 




I 
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I 













puter capable of reading words of up to 12 parallel bits. The mag- 
netometer words are broken into two, 12 bit words and time multiplexed 
on tape to satisfy the computer input conditions. Hence, for each 
sample, there are six 12 bit words to be processed. Even with the 
time speed up in reproduction the CDC 160 is operating below its 
capabilities. In an effort to utilize the l60 time, the data was 
put on tape, so that when it occured, the l60 would read at its maxi- 
mum rate and use the ”no data time'* to continue with preliiainary 
processing. Channels 13 and lU on the tape are input ready and dis- 
connect signals used to control the l60. And of course, for reason- 
able machine time efficiency, there is on hand enough tape to allow 
continuous recording for one week. 

The Ampex CPIOO was not designed specifically for this type of 
operation! hence the need for an interface to mate it with the 160. 

The interface must provide the following functions g 

a. Convert the pulse type data to D,C. level signals. 

b. Put the data on line when a request is supplied from the 
computer. 

c. Effectively disconnect itself from the lines when other 
equipment is called. 

d. Provide control information, i.e. input ready and discon- 
nect for the l60. 

e. Effectively isolate the computer from spurious transients 
from the tape machine. These normally occur dui'ing start - stop 
operations, but the recorder is somewhat susceptible to power line 
transients and they sometimes appear in its output. 
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Figures 8, 9 and 10 shows the block diagram of the interface 
provided* Fig* 11 is a photo of the preliminary processirig setup. 
The cards were assembled and mounted as a component on a laboratory 
cart, the tape unit is placed on the top of the cart and computer 
cable connectors provided. This eliminates the need for a permanent 
rack in the computer laboratoiy. 

The preliminary processing consists of the followings 

a, Reading the raw data, 

b. Checking the data and keeping a record of the errors* 
c* Writing the storage tape* 

d* Providing identifying blocks on the storage tape. 

e. Searching the storage tape for ary particular record or 
the end of the data so that one storage tape may be used for many 
week’s records, 

f. Direct display of the raw data* 

The data storage format used is illustrated in Fig* 7* 



STORAGE TAPE FORMAT 
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Figure 7 

With the time base used there are nearly 380li0 twelve bit words 
per 6,5 hours record. The 6*5 hours record is one 3600 roll of in- 
strumentation tape at 1 7/8 ips. Since the operation is stopped to 
change tape, it is convenient to put the identifying block at the 
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end of this irecord. Also, this number of data words can be process 
sed in the order of minutes by the CDCI60I4 so that if the computer 
is needed for other jobs| it can be stopped at the identifying blocks 
and return later to the same point with the proper index setting* 

The block length was set somewhat arbitrarily* The memory storage 
limit of the I60, write time of the I63 tape unit and suitable re= 
cords for handling in the l60i;, are the factors to be considered when 
chosing the block length. The instrumentation recorder cannot con-= 
veniently be stopped when writing the storage tapej hence one data 
point is lost when transferring data* This represents 22 points 
per 6.5 hours record or less than a minute in real time. The block 
length was chosen to be 3li56]_Q twelve bit xrords per block, the fig=- 
ure being corapatable xdth the above conditions. 

The identifying record consists of xrriting the month, day and 
year in flex code. The remaining cells of the block are loaded vrith 
zeros. It is then easy to separate the data from the identifying 
blocks since none of the data words will be all zeros. 

The fact t’lat the data had to be reconstructed in tl is Ic'Oi, x-jas 
also considered in chosing the block length. The information is 
buffered into the l60I^ in U8 bit words. The blocks consist of an 
even number of data points, i.e, six 12 bit or two li8 bit words* 

After a block is read into the machine and loaded into the A register, 
the magnetometer readings can be stored sequentially in memory with 
a loop of a few machine language instructions. The data analysis 
is then done using Fortran, alloxri.ng easy and versatile manipxalation 
of the data. The I60 and I60U input programs are included in Appcn=* 
dices II and III. 
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Reassembly is illustrated below 



160U WORDS 



Mag 1 (2U bits) 


Mag 2 (2U bits) j 




Mag 3 (2U bits) 


Mag 1 (2ii bits) 




Mag 2 (2U bits) 


Mag 3 (2U bits) 



Briefly, the operation is as follows? Magnetometer readings 
are recorded for long periods of time at a remote site on a protable 
magnetic tape units The tape unit is then returned to the process^ 
ing site, and using a time speed up of X 32, the intermediate pro- 
cessing is performed and the data stored* 

The storage tape along with suitable programs is then placed 
on file for processing by the computer facility at their convenience. 
The tape unit is then returned to magnetometer site for more records. 
The new data obtained is added to the storage tape so as not to de- 
stroy the existing data. 
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3. Data Analysis 



The experimental measurement program was set up to study the 
coherence between the magnetometersj i.e, the comparison of simulta- 
neous samples of the field, with the sensors located at different 
spatial locations. The coherence was to be examined with various 
local perturbations in the field. With the target detection objective 
in mind, the least complicated method of obtaining the target signal 
will result in the best instrumentation for field use. The problem 
is somewhat analogous to the radar problem where decisions must be 
made in the presence of noise, with associated probabilities of error 
in detection. 

There are several types of geomagnetic noise that occur. The 
most dominant is the diurnal variations. Superimposed on this are 
smaller variations known as micropulsations and noise caused by lo- 
cal disturbances. Sudden eruptions on the sun sometimes cause large 
fluctuations in the field. Several studies of the various noises 
have been made (3), (Ij.) and (^), For the proton precession magneto- 
meter, the instrument produces noise of major contribution to the 
spectra of small geomagnetic pulsations. For target detection, the 
signal of interest will most likely be of the same order of magni- 
tude as the small micropulsations and local distrubances. For geo- 
physical prospecting and mapping, the signals will be much larger. 

As with any raw data, there will be some points that are obvi- 
ously in error. The cause of errors is ustially uncertain, but some 
provision must be made in the processing to eliminate these errors 
so as not to distort the data. Care must be taken when performing 
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these operations so the target signals are not masked. Records of 
the errors should be kept in order that decisions about the useful- 
ness of the set of data can be made. 

The problem of coherence between the spatial samples requires 
that the instnament noise be determined^ with an ideal signal^ then 
measurements be conducted with no local perturbations in the field, 
and the samples compared. The simplest comparison is to plot the 
difference signal vs. time. With no local perturbation, this means 
essentially a plot of the system noise vs. time| hence the need for 
filtering is obvious. To use this method for target detection will 
require serious consideration of the filter time constants. 

The types of filtering 

1. A clipping filter (on the raw data) that clipped at ^ 200 y. 
The number of points outside of this range are recorded. This range 
was about twice the expected diurnal variation. The filter removes 
the errors from lost high order bits and other instrument soxu’ces. 

The typical error rat-e was about 1/50000. 

2. A simple three point average for display of the raw data on 
the DD 65 display xinit was used. This gives an operator a fast vis- 
ual display of the raw data. Decisions can then be made about the 
quality of the data. 

3. Eighteen point averaging, i.e. 36 seconds real time, gave 
partially smoothed 6,5 hoxirs graph records, 

U. A simple single section low pass filter was used of the 
form l/(s - a). The time constant was set at 100 sec, and 20 min, 
in real time. The implementation was done using 2 transform (6). 
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The filter design was done by Professor Harold A« Titus 



We have 



x„(s) 

Jln(S) 



a 



s - a 



which has a Z transform 

X(^(.2) / - e d 

In block diagram form this is seen in Fig« 29 




Figure 29 



In terms of the discrete time samples we have 

Xo(A-h)= /S ^ ^ ^vi ( 






— ccrvy^ 



T g sample period 
a s time constant 

The function was truxjkated at n § 20 for values of k ^ 20« 
The contribution to the filtered value from the terms multiplied 
by e"^^ for n ^ 20 is negligible. 
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The coherency Is a statistical parameterj and raist be repre= 
sented as such. 

Let? 



The differences of the unfiltered values were grouped in a frequency 
distribution which approximates the continuous density function for 
large sample size. Since the process is digital^ the allowable dif- 
ferences, in bits ares 0, ^ 1, ^ 2, ^ »«<■. Therefore, the frequency 
distribution is not a smooth cvirve. Since the instrument noise is 
randomly distributed, the envelope of the distribution is a normal 
density function. The mean was removed before plotting, and the 
statistic representing the standard deviation was computed according 



Comparison of the density function curve with an ideal signal 
input to that of magnetometer input gives a measure of the noise 
produced by the device. 

To f’urther define the coherence, the auto and cross correla- 
tion functions were computed and the coherency computed as defined s 



X = «1 - «2 

A 2 - ”1 - "3 
A 3 = »2 - »3 




where s N = the ntuaber of samples in a 6,5 hour record* 
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where s 



and 





L 



L 







This discrete computational form is 



and 




where s 

N g number of data values 
n 3 number of lags 
i - n 



Since the functions are from finite time series can only be made 

about 0.1 IJ. Tho 6.5 hour records will vary in form deperdirsg upon 
the time of day the records were made. Therefore, the correlation 



functions will not be the same shape for each reeord| however the 
coherence function as defined will be the samej and for perfectly 
coherent records is unity* The programs are included in Appendix 
III. 

It should be pointed out that when perturbations are placed in 
the field, the basic assumptions regarding correlation]; ioe* a sta^ 
tionary random process, are no longer applicable* 

Experimental data was taken to determine typical signatures for 
elementary perturbations. Fig® 12 shows the geometry of the detector 
centered coordinate system® 



PERTURBATION EXPERIMENT GEOMETH? 
1 - 




Figure 12 
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The dipole moment was orientated parallel to the earth's field and 
its motion was in the x direction such that dx/dt = constant# The 
approximation r 1 is valid# 

At the detectors 

H g ]j Ho^ = Hp^ 



Ands 



Hp 




where s 

M s dipole moment 

r s function of time and the geometry 
K s constant 

No attempt was made at this time to compare the measured values with 
the theoretical values ^ because the terrain does not lend itself 
keeping L, d, and v constant. The results are representative of 
typical signatures that X'Jould be found in practice. 
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U. Results 



l*.l Instrxanenta tloii a 

An ideal signal was substituted for the procession signal and 
set of records were made, as shown in Figs. 13 through The 
graph scales were left unchanged because it was desired to show the 
instrument errors on the same scale as the magnetic field scales. 

The magnetic field vs. time graphs were not included because the in- 
put signal was an audio oscillator and only a small drift puts it 
off the graph scale. The clock was 300 kc/s which makes 1 bit a 
0.36 / . 

Figure 28 shows the precesion signal from the various sensors. 
Figs. 26 and 27 are typical plots made from the analog section of the 
instrument. For comparison, the Rb®7 magnetometer output was adjust- 
ed to the same scale and plotted on the same graph. The Rb®^ magneto- 
meter output is an analog of frequency, while the proton magnetometer 
the analog is of period^ hence the reason for the opposite drift on 
the chart. Figure 27 is an analog record of the differences. The 
filter was a single section RC filter with time constant § 10 sec. 

On the miiLtiple curve plots showing coherence, the curves are identi- 
cal and cannot be distinguished from each other. These curves can 
be compared with the curves obtained from the sensor with no local 
pertTirbations in the field, for the contribution of the measuring 
instrument noise to the total noise. 

Difficulties with both instrumentation and computer tape have 
been encountered. The use of new computer tape for data storage and 
scratch tapes, and frequent head cleaning of the instrumentation 
recorder reduced the difficulties to a minimumc 
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Differences in magnetic intensity vs. 
time for an ideal signal input. 

DEL 1 s Hi - H2 
DEL 2 = Hi - H3 
DEL 3 r H2 - H3 

Filter a Z transform low pass with 
time constant of 20 minutes. 
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Normalized density function for the 
difference Hi - H3 - mean (Hi - H3)> ^ 

with an ideal signal input. ® 

Kean a 0.056 gamma 
Standard deviation r 0.72 



<s 

CS) 



I 




X-SCfl-E > s.08E+ea imtmnch. 
Y-SCfiLE - 2 ^-ai UNlT3/INCa 



flNDERSOM BOX 263 
DENSITY FUNCTION 



Figure I6 



X IN GAMMA Y 



IN FREQ 



33 




» >u T iHr-#.* ni K ncrruiyi TTtan^ 



Normalized density function for the 
difference Hp - H 3 - mean (H 2 - H 3 ), ^ 

with an ideal signal input. ® 

Kean 3 —0.0055 gamma 
Standard deviation = 0,72 
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Cross correlation for an ideal 
signal input. 

Record length ~ 6,5 hours 




m 



C31 635 0a<3 83;^ 



Figure 19 



^“SCfiLE - L6SE-f01 IMTS^IHCH 

y-scple - lesE+eii mira^iNCH. 

RNDERSON box 263 

croscorrelpition function y in products r in lrgs 



MO" zto 0Z0 csa SC0 f-jdo 



Coherence for an ideal signal input 
Record length s 6, S hours 
H]_2 = Coherence of Hx and H 2 
H13 = Coherence of Hi and H3 
H23 = Coherence of H2 and H3 
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Differences in magnetic intensity 
for an ideal signal input 
Record length = 1 hour 
Filter » 2 point average 
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X-SCALE = l.OOE + 00 UNITS/iNCH 
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FIRST DIFFERENCE IN MG FIELD T IN MINUTES H IN GAMMA 

Figure 21 
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Difference in magnetic intensity 
for an ideal signal input 
Record length s 1 hour 
Filter = Z transform low pass 
with time constant of 100 seconds. 
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Figure 22 
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Auto correlation for an ideal 
signal input. 

Record length - 1 hour 
H]_ r Correlation (Hi - mean Hi) 
H2 s Correlation (H2 - mean H2) 
H3 = Correlation (H2 - mean H3) 
Lag - A t 
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Figure 23 
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Cross correlation for an ideal signal input. 
Record length 3 1 hour 

Hi2 s Correlation (Hi - mean Hi and H2 - mean H2) 

H13 m Correlation (Hi - mean Hi and H3 - mean H3) 

H23 = Correlation (H2 - mean H2 and H3 - mean H3) 

Las a A t 
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Coherence for an ideal signal input 

Record length 3 1 hour 

Hi 2 - Coherence of Hi and H2 

H13 a Coherence of Hi and H 3 

H23 = Coherence of H 2 and H 3 
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TYPICAL ANALOG PLOT OF 
MAGNETIC INTENSITY VS. TIME 
FILTER . RC WITH TIME CONSTANT OF 6 SECONDS 
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TYPICAL ANALOG PLOT OF 
THE DIFFERENCES VS. TIME 
FILTER » RC WITH TIME CONSTANT OF 6 SECONDS 
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li«2 Magnetic field moasvtreiaents 

Only representative results of the magnetic field measuranients 
are presented^ since the measurements are not complete at this time* 
Samples of plots riade with the clock at 100 kc/s are included to 
illustrate the effect of the clock frequency on the results* For 
the associated density function of the differences, the allowable 
values are 0 ^ 1*1, ^ 2*2, • ® <> in gamma, v;ith the clock at 100 

kc/s. The expanded difference vs. time shows clearly the ^ last 
bit limitation* Comparing the density functions with the plots for 
the clock at 300 kc, shov/s that most of the noise associated wi.th 
the measurement is due to the last bit and can be reduced ty increas- 
ing the clock frequency* The perturbation experiments were conducted 
with the clock at 300 kc/s* The resuirs of the perturoauxon experi- 
ments show clearly the effect of the digital filtering, and the im- 
poirbance of rime constant consideration for target detection. The 
oscillator is to be run 1 mc/s for further experiments. Also, inapy 
of the errors that appear in the data were the direct result of the 
instrumentation tape. Mew tape has been purchased and should elim- 
inate most of these errors, in future records. 

As pointed out previously, the basic assumptions regarding cor- 
relation are not valid, vrhen the noise is not stationai-yj however, 
since correlation detectors can be instrumented and may prove useful 
in the detection problem. Some of the eon'elation function plots are 
included to illustrate the effect of local pertui'bations on the func- 
tions. 
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Conclusions 



For target detections where the signals under consideration are 
in the ^ gammas or less ranges it will be difficult to detect the 
target •without resorting to sophisticated signal processing* In 
this range of signalss consideration should be given to other instru= 
mentsj such as the Rubidium optically pumped luagne tome ter s whose 
theoretical sensitivi-ty is OoOOl gammas especially where the detector 
is stationary’s orientation is not a problen* It can be demonstrated 
that when using the same measurement technique as for the proton 
magnetometers the instrument and processing noise will be an order of 
magnitude less than the theoretical limits imposed by the sensiti- 
vity of the instrument. 

For signals in the J to 1^ gamuia ranges detection could be ac" 
complished with analog type correlation detectors. These could be 
implemen'bed with operational amplifiers. 

Signals above this ran^e could be detected with the simple dif- 
ferencing techniques and the associated RcC, filters. 

Increasing the clock frequency and averaging the period over 
more than 1G2U cycles would reduce the laeasm’ement noise to a smaller 
part of the to’tal noise. The target signals could then be reduced 
by a factor of three^ with no clsange in the Instrujr.en'bation. 

Full utilization of the inherently narrow line width of proton 
magnetometer cannot be achieved using the standard measuring teeh= 
nique. Other techniques have been suggested by Professor Carl I. 
Kennekenj such as utilising a rapid coherent polarization in a phas©= 
lock loop arrangement* Further devolopiuent work needs to be don^ 
to explore these possibilitieso 
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Typical earth’s magnetic field vs. 
time for clock at 100 kc/s and no 
local perturbations. 

Record length r 6,5 hours 
Hi s Hi - mean Hi 
H 2 s H 2 “ mean Hi 
H^ r H 3 - mean Hi 

Filter s 2 transform low pass with 
time constant a 20 minutes 
Kean Hi a 0.50917 E + 05 
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Difference in magnetic intensity vs. 
time for clock at 100 kc/s and no 
local perturbations. 

Record length s 6.5 hours 
DSL 1 3 Hi - H2 

DEL 2 r Hi - 

DEL 3 s H2 •• 

Filter a 13 point average 
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Dirrerence in magnetic intensity vs* 
time for clock at 100 kc/s and no 
local perturbations. 

Record length =6,5 hours 
DEL 1 . Hi - H2 
DEL 2 « Hi - H 3 
DSL 3 s H 2 - H 3 

Filter = 2 transform lov; pass with 
time change « 20 minutes 
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Normalized density function of the 
difference, D. No local perturbations 
Record length = 6,5 hours 
Clock - 100 kc/s ® 

(Hi - H 2 ) - mean (Hi - H 2 ) - d 
' Mean m -3.21 gamma 
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Normalized density function of the 

difference, D» No local pertuii>ations 

Record length =6,5 hours ^ 

Clock s 100 kc/s 

(Hi - H 3 ) - mean (Hi - H 3 ) = D 

Mean a 3,16 gamma 
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Auto correlation for earth's magnetic 

field vs. time. No local pertuit>ations 

Record length s 6.5 hours 

Clock s 100 kc/s 

H]_ s Correlation (Hi - mean Hi) 

H2 = Correlation (H2 - mean H2) 

H3 s Correlation (H3 - mean H3) 

Lag * A t 
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Cross correlation for earth’s magnetic 
field vs. time. No local perturbations 
Record length a 6,5 hours 
Clock = 100 kc/s 

Hi2 r Correlation (Hi - mean PJi and H2 mean H2) 

H13 - Correlation (Hi - mean Hi and H3 - mean H3) 

H23 5 Correlation (H2 - mean H2 and H3 - mean H3) 

Lag -At 
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Coherence for earth's magnetic field 
vso time. No local perturbations 
Record length r 6.5 hours 
Clock s 100 kc/s 

“ Coherence of and H 2 
^13 = Coherence of Hi and H 3 
H 23 = Coherence of H 2 and H 3 
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Normalized density function of the 
differencCjD. N3 local perturbations 
Record length r 6,5 hoxirs 
Clock s 300 kc/s 

(Hi - H2) - mean (Hi - H2) = D 
Mean = -li«09 gamma 
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Nornialized density function of the 
differencejD. local perturbations 
Record length =6,5 hours 
Clock = 300 kc/s 

(^2 ” ^ 3 ) ~ mean (H 2 ~ H 3 ) - d 
M ean - -0,33 gamma 
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Cross correlation for earth’s magnetic 
field vs. time. No local perturbations 
Record length s 6.5 hours 
Clock s 300 kc/s 

Hi 2 2 Correlation (Hx - mean Hi and - 

H13 = Correlation (Hi - mean and H3 - 

H23 - Correlation (H2 - mean H2 and - 

Lag s A t 



mean H2) 
mean H3) 
mean H3) 
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Coherence for earth’s magnetic field 

vs. time. No local perturbations 

Record length s 6.5 hours 

Clock 3 300 kc/s 

Hi 2 = Coherence of and H2 

H13 = Coherence of and H3 

H23 = Coherence of H2 and H3 
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Earth's magnetic field vs. time with 
20 gamma step function applied to 
Record length r 6.5 hours 
Clock r 300 kc/s 

Filter = Z transform low pass with 

time constant of 20 minutes 

Mean = (mean + mean H2 -V mean H3)/3 

% s Ph - Mean 

H2 = K2 - Mean 

H3 : H, - Mean 

Mean r 0.509307 E + 05 
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Earth's magnetic field vs, time 
Filter - 18 point average 
Clock r 300 kc/s 
Record length a 6.5 hours 
Hx = H]_ - mean Hx 
H2 s H2 - msan Hx 
H 3 = H 3 - mean Hx 
Perturoations : 

5 gamma and +• 10 gamma applied 
at constant velocity to H 3 , - 2 
gamma, - 5 gairma, - 5 gamma applied 
at constant velocity to H 3 , H 2 and 
Hx during first hour, A 5 gamma 
step applied to H 3 at t a 1 hour and 
removed at t ^ 1 a/U hours. 
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Earth* s magnetic field vs, time 
Filter = 2 point average 
Clock 3 300 kc/s 
Record length = 1 hour 
% = % " mean 
^2 z Hp - mean H]_ 

H3 r Ho - mean H;j_ 

Per turoa tions : 

+ 5 gamma and ± 10 gamma applied 
at constant velocity to H3, -2 
gamma, -5 gamma, -5 gamma applied at 
constant velocity to Ho, H2 and H^ 
during first hour, A 5 gamma step 
applied to H3 at t = 1 hour and 
removed at 1 3/U hours. 
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EARTHS MAGKTJTIC FIELD VS TIME T IN MINUTES H IN GAK-IA 



I 

[ 
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Earth’s magnetic field vs. time 
Filter r Z transform low pass with 
time constant of 100 seconds. ' 
Clock s 300 kc/s 
Record length = 1 hour 

^1 - ^1 " Jaean 

^2 “ '^2 ” % 

H3 = H3 - mean 

Perturbations : 

- 5 gamm and - 10 gamma applied 
at constant velocity to Ho. - 2 
gamma, - 5 gamma, - 5 gamma applied 
at constant velocity to H3, H, and 
H], during first hour. A 5 gamma 
step applied to Ho at t - 1 hour 
and removed at t a 1 3A"*ho\irs. 




X SCALE - l.OOE-OO UNITS/iNCH 
I SCALE - l.OOE-Ol units/inch 

earths MAGNETIC FIELD VS TIME T IN MINUTES H IN GAM-IA 



1-32 




= Correlation 
H2 5 Correlation 



Auto correlation for earth's magnetic 
field vs. time. 

Record length a 1 hour 
Clock s 300 kc/s 

- mean Hi) 

(H2 - mean H, , 

- Correlation (K-? - mean nf) 

ligsAt ^ ^ 

Perturbations'; 

* 5 gamma and i 10 gamma applied at 
constant velocity to H,. -2 gamma, 

-5 gamma, -5 gamma applied at constant 
velocity to H^, H2 snd H]_ during first 
hour, A 5 gamma step applied to K3 at 
t = 1 hour and removed at t s 1 3/h hovirs. 



oai 



002 



003 



00^ 



005 



0oe 



K"5CflLE “ L00E+01 UNITS/INCH. 
Y -SCALE » 1.00E+00 UNITS/INCH. 



IDERSOM BOX 263 

ITOCORRELATION FUfiCTION Y IM PRODUCTS 



m 
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Cross correlation for earth's magnetic 
field vs. time. 

Record length - 1 hour 
Clock = 300 kc 7 s 
Hi 2 r Correlation 
Hx3 = Correlation 
H23 = Correlation 
Lag r A t- 
Perturbations; 

i 5 gamma and i 10 gamma applied at 
constant velocity to Ko. -2 gamma, 

-5 gamma, -5 gamma applied at constant 
velocity to H3, Hp and during first 
hour. A 5 gamma step applied to Ho at 
t a 1 hour and removed at t r 1 3 /U hours. 



H13 



(Ht -mean and Hp -mean H2) 
(Hf -mean and Ho -mean H^) 
(H2 -mean Hp and -mean H^) 



f 63l 002 803 805 ^5 006 O^ 



X -SCALE - LOOE+01 UNITS/INCH. 

Y-SC^LE - 1.00E+00 UNITS/INCH. 

3ERS0M BOX 263 

OSCORRELfiTIOH FUNCTION Y IN PRODUCTS X IN LAGS 
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m 



H. 






Coherence for earth’s magnetic 
vs « time . 

Record length ~ 1 hour 
Clock - 300 kc/s 
Ht 2 : Coherence of and 

^13 = Coherence of and 

H23 r Coherence of H2 and Ho 

Lag r At 
Perturbations ; 

± S gamma and ± 10 gamma applied at 
constant velocity to Ho. -2 gamma, 

-5 gamma, -5 gamma applied at constant 
velocity to H 3 , Ho and H;i_ during first 
•hour. A 5 gamma step applied to at 
t r 1 hour and removed at t ^ L 3 /u hours. 




881 



882 



883 



50’i 



oas 



08Q 



08<7 



- 1.03E+01 IMTS^IMCH. 
Y-5CflLE » 2.08E-81 UNITS/IHi^ 



ANDERSON BOX 263 

COHERENCE FUNCTION Y IN COH X IN LAGS 
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160 FTOGRAMS AllD 
OPERATOR IHSTRUCTIOI'iS 



II 



mmS 















! 



3. 
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160 PROGRAMS AND OPERATOR INSTRUCTIONS 



Starting; Addresses 



Function 



0000 

0002 

0004 

0006 

0010 

0100 



Rewind 

Read 

Search 

Write 

Write identification Block 
Main Input Program 



Error Stops 



Cause 



ERROl 

ERR02 

ERR04 

ERR03 



773 of current input in error 
Memory storage not complete 
End of storage tape 
163 not ready 



Addresses of 
Important Counters 



Information 



0057 

0045 

0046 



Number of blocks v/ritten on 
storage tape 

Current storage location in 
memory 

Number of zero words in mem- 
ory 



Operator Instructions 

1. Load program and make the necessary cable connections o 

2. Set Identifying records in cells 3O3 to 370. 

3. Write identifying block (run from P s 010) « Stop at 
P s 100. 

4. Set cell 100 s 2200 

57 s 0 

5. Start Ampex C«P, 100 and run from P g 100, Stop when 
tape is finished or ERR stop. 
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0100 

0101 

0102 

0105 

0104 

0105 

0106 

0107 

0110 

0111 

0112 

0115 

0114 

0115 

0116 

0117 

0120 

0121 

0122 

0125 

0124 

0125 

0126 

0127 

0150 

0151 

0152 

0155 

0154 

0155 

0156 

0157 

0140 

0141 

0142 

0145 

0144 

0145 

0146 

0147 

0150 

0151 

0152 

0153 

0154 

0155 

0156 
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2200 


LDFOO 




1000 


1000 




4045 


STD45 


(ADDR) 


0400 


LDNOO 


4046 


STD46 


(GNTRl) 


0400 


LDNOO 


4150 


ST 1 50 




4151 


ST 1 51 




4152 


ST 1 52 




4155 


ST 1 53 




4154 


ST 1 54 




4155 


ST 1 55 




7500 


SXFOO 




0500 


0500 




7205 


INPO 3 




0620 


0620 




6102 


NZF02 




0600 


0600 




1600 


LSFOO 




0607 


ADNO 7 




6115 


NZFI 5 




2150 


LDI 50 




6015 


ZJFI 5 




2151 


LDI 5 I 




6011 


ZJFll 




2152 


LDI 52 




6007 


ZJF 07 




2153 


LDI 55 




6005 


ZJFO 5 




2154 


LDI 54 




6005 


ZJFO 3 




2155 


LDI 55 




6105 


NZFO 5 




5446 


A0D46 




077Y 


SBN77 




6536 


NZB36 




0001 


EElROl 




2050 


LBD50 




4047 


STD47 




2147 


LDI47 




4145 


ST 145 




5445 


A0D45 




5447 


A0D47 




3456 


SBD56 




6505 


NZB 05 




2200 


LDFOO 




4777 


4777 





II=.5 



0157 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0170 

0171 

0044 

0045 

0046 

0047 

0050 

0051 

0052 

0053 

0054 

0055 

0057 
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3445 

6653 

2044 

3445 

6003 

6657 

0002 

7101 

0172 

1000 

7700 



SBD45 

PJB53 

LDD44 

SBD45 

ZJF03 

PJB57 

EEIR02 

JFIOl 

WRITE 



HLTOO 



ADDITIONAL CELLS USED 



CNTR 1 
CNTR 2 

BUFFER STORAGE 
BUFFER STORAGE 
BUFFER STORAGE 
BUFFER STORAGE 
BUFFER STORAGE 
BUFFER STORAGE 

NUI4BER OF MEMORY DUMPS 



LIMIT 

ADDR 



(PUT IN LIMIT) 
7600 
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0172 

0173 

0174 

0175 

0176 

0177 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0210 

0211 

0212 

0213 

0214 

0215 

0216 

0217 

0220 

0221 

0222 

0223 

0163 

0164 

0165 

0166 

0167 

0170 

0171 
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5457 


A0D57 


2600 


LCFOO 


7776 

4060 


STD60 


7500 


EXFOO 


1141 


REQ STAT 


7600 


INAOO 


0202 


LPN02 


6004 


ZJF04 


5460 


A0D60 


6506 


NZB06 


0003 


ERR03 


7500 


EXFOO 


2111 


2111 


7303 


0UT03 


7600 


TEEIM OUT (LI 


6102 


NZF02 


1000 


INIT WRIT 


7500 


EXFOO 


1141 


REQ STAT 


7600 


INAOO 


0240 


LPN40 


6103 


NZF03 


7101 


JFIOl 


0100 


00 


7700 


HLTOO 



+ 1 ) 



DUI'^P TO PAPER TAPE 



7500 

4104 

7303 

7600 

6102 

1000 

7700 
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DAILY IDENTIFICATION 



USE 8 CELLS 
0030 


0070 


EXAMPLE 

3 


0031 


0064 


0032 


0032 


J 


0033 


0030 


A month 


0034 


0060 


N 


0035 


0072 


6 


0036 


0062 


4 year 


0037 


0045 


CAR RET 



11=7 
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S'p/JRT’ 



DAILY RECORD TO M.T. 



Up£> 
S O 



S'PX 

looo 




0224 

0225 

0226 

0227 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0260 

0261 

0262 

0263 

0264 

0265 

0266 

0267 

0270 

0271 

0272 

0273 

0274 
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DAILY RECORD TO MoT, 



2030 


LDD30 


4100 


ST 100 


1000 


1000 


2031 


LDD31 


4100 


ST 100 


1001 


1001 


2032 


LDD32 


4100 


ST 100 


1002 


1002 


2033 


LDD33 


4100 


STIOO 


1003 


1003 


2034 


LDD34 


4100 


STIOO 


1004 


1004 


2035 


LDD35 


4100 


STIOO 


1005 


1005 


2036 


LDD36 


4100 


STIOO 


1006 


1006 


2037 


LDD37 


4100 


STIOO 


1007 


1007 


2200 


LDFOO 


7700 


7700 


4100 


STIOO 


0100 


0100 


2200 


LDFOO 


1010 


1010 


4045 


STD45 (ADDR) 


0400 


LDNOO 


4145 


ST 1 45 


2044 


LDD44 


3445 


SBD45 


6003 


ZJF03 


5445 


A0D45 (ADDR) 


6506 


NZBO6 


7101 


JFIOl 


0172 


0172 (DU14P) 


7700 


NORl-IAL STOP 
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REVJIND, READ AND 







rki i. jL^ 



Rmm). RFAD AND 
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Start 



Rewind 0275 


0400 


LDNOO 


0276 


4043 


STD43 


0277 


7500 


EXFOO 


0300 


1141 


REQ STAT 


0301 


7600 


INAOO 


0302 


0202 


LPN02 


0303 


6504 


NZB04 


0304 


7500 


EXFOO 


0305 


1161 


REWIND 


0306 


7500 


EXFOO 


0307 


1141 


REQ STAT 


0310 


7600 


INAOO 


0311 


0202 


LPN02 


0312 


6504 


NZB04 


Read 0313 


7700 


HLTOO 


Start 0314 


7500 


EXFOO 


0315 


1141 


REQ STAT 


0316 


7600 


INAOC 


0317 


0202 


LPN02 


0320 


6504 


NZB04 


0321 


7500 


EXFOO 


0322 


2111 


EXFll READ 


0323 


7203 


0UT03 


0324 


7600 


TERM READ 


0325 


6102 


NZF02 


0326 


1000 


INIT READ 


0327 


7500 


EXFOO 


0330 


1141 


REQ STAT 


0331 


7600 


INAOO 


0332 


0240 


LPN40 


0333 


6002 


ZJF02 


0334 


0004 


ERR04 (END 


0335 


2043 


LDD43 TA; 


0336 


0701 


SBNOl 


0337 


6004 


ZJF04 


0340 


6525 


NZB25 


Search 0341 


0401 


LDNOl 


Start 0342 


4043 


STD43 ( SEA] 


0343 


2100 


LDIOO 


0344 


1000 


1000 


0345 


3430 


SBD30 


0346 


6532 


NZB32 


0347 


2100 


LDIOO 


0350 


1001 


1001 


0351 


3431 


SBD31 


0352 


6536 


NZB36 


0353 


2100 


LDIOO 
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0354 


1002 






1002 


0355 


3432 






SBD32 


0356 


6542 






NZB42 


0357 


2100 






LDIOO 


0360 


1003 






1003 


0361 


3433 






SBD33 


0362 


6546 






NZB46 


0363 


2100 






LDIOO 


0364 


1004 






1004 


0365 


3434 






SBD34 


0366 


6552 






N2B52 


0367 


2100 






LDIOO 


0370 


1005 






1005 


0371 


3435 






SBD35 


0372 


6556 






NZB56 


0373 


2100 






LDIOO 


0374 


1006 






1006 


0375 


3436 






SBD36 


0376 


6114 






NZF14 


0377 


2100 






LDIOO 


0400 


1007 






1007 


0401 


3437 






SBD37 


0402 


6110 






NZFIO 


0403 


2100 






LDIOO 


0404 


1010 






1010 


0405 


6105 






NZFO 5 


0406 


0400 






LDNOO 


0407 


4043 






STD43 


0410 


7101 






JFIOl 


0411 


0077 






0077 


0412 


7101 






JFIOl 


0413 


0314 










TEST TAPE 22 BLOCKS 


DATA 




M * 4714510 


LIMIT 


CELL 


44 = 


7600 




ADDR 


CELL 


45 




» 0013 4051r 


CNTR 2 


CELL 


21 




N2 Ml 


Ml 


CELL 


22 






N2 


CELL 


23 






CNTR 1 


CELL 


20 




0526 


7700 






HLTOO 


0527 


0400 






LDNOO 


0530 


4021 






STD21 


0531 


2200 






LDFOO 


0532 


1000 






1000 
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0533 

053^ 

0535 

0536 

0537 

0540 

0541 

0542 

0543 

0544 

0545 

0546 

0547 

0550 

0551 

0552 

0553 

0554 

0555 

0556 

0557 

0560 

0561 

0562 

0563 

0564 

0565 

0566 

0567 

0570 

0571 

0572 

0573 
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4045 

0413 

4023 

0400 

4020 

2200 

4051 

4022 

2021 

0726 

6417 

2023 

4145 

5445 

2022 

4145 

5445 

5422 

5420 
2200 
4777 
3445 

6207 

2044 

3445 

6104 

5421 
7101 
0172 
2020 
0777 
6434 

6525 



STD45 

LDN13 

STD23 

LDNOO 

STD20 

LDFOO 

4051 
STD22 
LDD21 
PPN26 
Z0B17 
LDD23 
STI45 
A0D45 
LDD22 
ST 1 45 
A0D45 
A0D22 
A0D20 
LDFOO 

4777 

SBD45 

PJFO7 

LDD44 

SBD45 

NZF04 

A0D21 

JFIOl 

0172 

LDD20 

SBN77 

ZJB34 

NZB25 
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WRITE TEST TAPE 
22 BLOCKS DATA 















f 

I 
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FORTRAN PROGRAMS 
ANALYSIS PROGRIIi 
SATELLITE DISPUY PRXRAM 



UIMtNSlCJN NDATA (664) , DATE( 00 ), H( 17 28) , HONE ( 900 , IT I TLE ( 1 2) , 
1HTW0(900) ,HTHRE (900) , T ( 9CC 1 , 8 IN 1 ( 10 0 ) » E I N2 ( 1 00 ) , B I N3 ( 1 00 ) , X ( 1 00 ) 
no 5 J = 1,14 
READ 5,DATi-(J) 

5 E0RNAT(C16) 

COr>!(NA.SK = 000000007777777 j'B) 

ENA(KDATA) , INA( 1 ) . 

STA( INI D , INA(664) . 

STA ( ITlRN) ,LNI ( 0) . 

6 CALL READ(:>, INI T, ITERK, 1 ) 

DO 70 1=3,.- 

SLSl (7) 

7 AOATE = DATE ( I ) 

BDATE = DATE ( 1+ 1 ) 

LDA( ADATE) ,SUR( NDATA +1 ) , AJPl (6) , 

LDA(BDATE) ,SU8( NDATA +2), AJPl (6). 

NBLOC = 0 

30 CALL READ(o, INI T, ITERM, 1 ) 

NZRC = 0 
DO 81 J=1 ,864 
IF(NDATA(J)) 82,82,81 
82 NZRO = NZRO+1 
81 CONTINUE 

NO = NDATA! 5) 

LDA(ND) ,AJP0( 10 ). 

K = C 

DO 40 J=1 ,864 
K = K+1 

INDATA = NDATA! J) 

LDA( INDATA) , LDG ( MASK ) , ARS { 4' 4 ) ,STL(INLIST). 

H(K)=INLIST 
K = K+1 

LDA( INDATA) , LDQ (MASK ),STL(1NLIST). 

SLSl (40 . 

40 H(K)=INL1ST 
AVt=C. 

DO 50 J =1 , 1728 
H(J) = H(J)/3. 

H( J ) = 1024./H( J ) *2348400 .0 
50 AVE = H(J) + AVE 



PRINT 26, AVE 
C ZAPP FILTER 

NERR = 0 

AVE = 0. 

HBAR = 50950.0 
DO 51 J=1 ,1728 
DEV = A8SF(H( J) -HBAR) 

IFIDEV-200.) 51,52,52 
JCGUNT = J 

52 NERR = NERR + 1 
IFIJ-3) 53,53,54 

53 H( J) = H( J+3) 

GO TO 51 

54 H(J) = H(J-3) 

51 AVE = AVE + H(J) 

AVE = AVE/1728. 

PRINT 26, AVE 
PRINT 97, NERR, NZRO, JCOUNT 
NBLOC = NBLOC + 1 
WRITE TAPE 2, H 
GO TO 30 
10 REWIND 2 

C THIS SECTION COMPUTES H VS TIME 
C 18 PT AVE N IS TOTAL NUMBER CF POINTS PER DAY 
TIME = - 6.5/704. 

DO 60 L=l, NBLOC 
READ TAPE 2,H 
DO 62 J=1 ,32 
HIBAR = 0. 

H2BAR = 0. 

H3RAR = 0. 

DO 61 K=1 ,16 

M = 3*K - 2 + ( J - 1 )*5 4 III-l 

HIBAR = H(M) + HIBAR 
H2BAR = H(M+1) + H2BAR 

61 H3BAR = H(M+2) + H3BAR 

N = J + (L-1 )*32 

HONE(N) = Hi BAR/18. 

HTWC(N) = H2 BAR/ 1-8. 

HTHRE(N) = H3BAR/18. 

TIME = TIME + 6.5/704. 

62 T(N) = TIME 

r 



RE^^IND 2 
PR 1. N'T 2 7, N 
PRINT 27, M 
NPTS=N 

C MEAN VALUE OF DIFFERENCES 
JERR = 0 

I ERR = 0 APPENDIX III 

NERR = 0 

DELI = 0. 

0EL2 = 0. 

DEL3 = 0. 

DO 64 J = 1,NBLCC 
READ TAPE 2,h 
DO 6U L = 1,576 
M = 3*L-2 

DEL = H(M + 1) - HIM) 

IF! ABSF (DEL)- UO. ) 65,66,66 

65 DELI = DELI + DEL 
GO TO 67 

66 NERR = NERR + 1 

67 DEL = HIM + 1 ) -HIM + 2 ) 

IF(ABSFIDEL) - 40.) 68,69,69 

68 0EL2 = DEL2 + DEL 
GO TO 71 

69 I ERR = lERR 4 1 

71 DEL = HIM) - HIM + 2) 

IFIABSFIDEL) - 40.) 72,73,73 

72 DEL3 = DEL:i 4 DEL 
GO TO 64 

73 JERR = JERR 4 1 
64 CONTINUE 

REWIND 2 

PRINT 97, NLRR, I ERR, JERR 
97 FORMAT! 3120 
ABLGC = NBLOC 
DELI = DELI / (576. *AELOC ) 

DEL2 = DEL2/ 1576. «A8L0C ) 

DEL3 = DELS/ 1576. »ABLOC ) 

PRINT 49, DELI ,DEL2, DELS 
C DENSITY FUNCTION 
STDEV1=0. 

ST0EV2=0. 

STDEV3=0. 

NERR = 0 
lERR = 0 
JERR = 0 
XL=-25. 

DO 45 K=l,100 
B I N 1 I K ) =C . 

BIN2(K)=0. 

BIN3I K) =0. 

X(K)=XL 
45 XL= XL40.5 

DO 74 J=l, NBLOC 
READ TAPE 2,H 
DO 74 L = 1,576 
M = 3»L-2 

STDEV1= IHIM41 )-HIM)-CEL1 )*o2 4STDEV1 
STDEV2= (H(M4l)-HIM42)-CEL )**2 4 STDEV2 
STU£V3= (H(M)-HIM42)-DEL5)»*2 4STDEV3 
DEL = HIM41 ) - HIM) 

IF! ABSF I DEL-DEL D-25. ) 7 5, 75, 76 

76 NERR = NERR 4 1 

75 DEL = HIM41 ) -HIM42) 

IF! ABSF I DEL-DEL 2) -25. ) 7 7, 77, 78 
78 IFRR = lERR 4 1 

77 DEL = HIM) - HIM42) 

IF! ABSF I DEL-DEL 3) -25. ) 2 9, 29,80 
80 JERR = JERR 4 1 
29 DO 31 K=1 ,3 

IF(K-2)32,33,34 

32 DEL=H(M4l )-H(M) -DELI 
GO TO 35 

33 DEL=HIM4l )-H(M42)-DEL2 
GO TO 35 

34 DEL=H(M)-HIM42)-DEL3 

35 XL= -25. 

DO 36 N=1 , 1 00 

IFIDEL-XL) 37,37,38 III-2 

38 IF(DEL-(XL40.5) )39,39,37 

39 IFlK-2) 42,43,44 

42 PINl (N)=BIN1 (N)4l . 

GO TO 31 

43 D:N2(N) = BIN2(N)4l. 

GO TO 31 
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O O U U K I 1 l\ u t 
31 CCNTINbE 
7U COMINUh 

STbEVl =SQRTF(STDEV1/{576.=ABL0C) ) 

STDEV2 =SQRTF(ST0EV2/(576.«ARL0C) ) 

STDEV3 =SQRTF(STDEV5/ ( 57d.« A6L0C) ) 

PRINT Ug.STDEV 1 ,STDEV2, STDfcV3 
PRINT 97,NERR, lERR, JERR 

remind 2 

NORMALIZED DENSITY FUNCTION 
AMAX=BIN1 ( 1 ) 

BMAX=BIN2{ 1 ) 

CMAX=RIN3{ 1 ) 

no 200 J=1 , 1 00 

IF ( AM AX-BIN 1 (J) )201 ,202,20 

201 AMAX=RIN1{J) 

202 IF! BMAX-eiN2( J) )203,204,20L 

(-) + + t + * + J + + + t + S + t ♦ ++=!= + //++ + +++ + + + 

20U IF (CMAX-BIN3 { J) 1205,200,200 
205 CMAX=BIN3(J) 

200 CONTI NUN 



// 



PRINT U9, AMAX, BMAX,CMAX 
DO 206 J=1, 100 
BIN1 ( J) =BIN1 ( J) /AMAX 
BIN2( J)=DIN2( JI/BMAX 
206 BIN3( J)=BIN3( J) /CMAX 
DO 46 K = 1 , 12 
46 ITITLE(K)=8H 

ITITLEI 1 )=CHANO£RSCN 
ITITLE(2)=8H BOX 2o3 
ITITLE{7)=£H DENSITY 
ITITL£(8)=aH FUNCTIO 
IT ITLE(9)=6HN X in G 
ITITLEI 10)=8HAMMA Y 
ITI TLEII 1 )=8HIN FREC 
LA6A=4H pi 
LAR 8=4H P2 
LAbC=4H P3 

CALL DRAW(100,X,BIN1 ,0,C,LAbA, ITI TLE, 0,0, 0,4, 0,2, 8, 8,0, LA ST) 
CALL 0RAW(|00,X,8IN2,0,C,LABB, I T ITL E , 0 , 0 , 0 , 4 , 0 , 2 , 8 , 8 , 0 , LA S T ) 
CALL DRAW! 100,X,BIN3,0, C,LABC, I TITLE, 0,0, 0,4, 0,2, 8, 8,0, LA ST) 
506 N = N PTS 

H0NEAV=0. 

HTWCAV=0. 



HTHREAV=0. 

DO 48 L=1,N 
HONEAV=HONEAV+HCNE{ L ) 
HTWCAV=HTWf;AV+H.TWC(L ) 

48 HTHREAV=HTHREAV+HTHRE(L ) 

SV04 = N 

HONEAV = H0NEAV/SV04 

HTWCAV = HTWOAV/SVOU 

HTHREAV = HTHREAV/SV04 

PRINT 49, HONEAV, HTWOAV, HTHREAV 

49 FORMAT! 3E20. 7) 

PRINT 27,NDLCC 

27 FORMAT! 17) 

HONEAV= IHONEAV +HTWOAV +HTHREAV1/3. 
PRINT 26, HONEAV 
26 FORMAT! E20. 7) 

WRITE TAPE 2, HONE 
WRITE TAPE 2, HTWO 
WRITE TAPE 2, HTHRE 
REWIND 2 
DO 90 L = 1,N 

HONEIL) = HONEIL) - HONEAV 
HTWCIL) = HTWOIL) - HONEAV 
90 HTHREIL) = HTHREIL) - HCNEAV 
DO 1 L=1,12 
1 ITITLEIL) = 8H 

ITITLEI 1 )=8HANDERS0N 
ITITLE!2)=8H BOX 263 
ITITLE!7)=3H EARTHS 
ITITLE!8)=8HMAGNETIC 
ITITLE!9)=SH FIELD V 
IT ITLE! 10) = 8HS TIME 
ITITLEI 1 1 )=8HT IN HRS 

ITITLEI 12)=8H H GAMMA III-3 

LABA=4H HI 
LABB=4H H2 
LABC=4H H3 



C Z 




transform 

1. I.ILIA/inn 



LAST ) 
LAST ) 
LAST ) 



I 




i 

.1 



IF(J-20) 302,303,303 
502 M=J 

GO TO soil 

303 M=;o 

304 00 301 K=1,M APPENDIX III 

L= K-1 



DL=L 

ASUI^ = E:XPF( -AI *PL 1+ASUP 
BSUM =1: XPF(- AT *6L) IP HONE (J-D + BSbK 
CSUK =EXPF(-AT*eL) «HTWO(J-L ) + CSUM 
301 OSUP = EXPr ( -A1 «DL) iiHTHRE( J-D+OSU'-l 
HONE( J)=BSUM/ASUFi 
HTWO(J)= CSUN'/ASUM 
300 HTHRE(J1= DSUM/ASUM 
ITITLE(2)=BH FILTER 

CALL DRAW(N, T,HCNE, 1,0,LABA,ITITL£, l.C,10.0,5,1,2,2,8,1C,0,LAST) 
CALL ORAW(N, T,HTW0,2, 0, LABB,ITITLE, 1.0, 10. 0,5, 1,2, 2, 8, 10,0, LA ST) 
CALL DRAW (N, T,HTHRE, 3,0 ,LABC, ITITLE, 1.0, 10., 5, 1,2, 2, 8, 10,0, LAST) 
C DIFFERENCES OF FILTERED VALVES 
DO 306 J = 1 ,N 
H( 4)= HONE( J ) 

HONE( J)= H( J)- HTWO( J) 

STORE = HTWO(J) 

HTWO(J)= H( J )-HTHRE( J ) 

306 HTHRL(J)= STCRE- HTHRE(J) 

ITITLE(7)=6HFILTERED 
ITITLEt 8)=eH0IFFEREN 
TTITLE(9)=8HCE IN NA 
ITITLEI 10)=8HG FIELD 
LABA=4HD£L1 
LABB=4HDEL2 



LABC=4HDEL3 

CALL DRAW(N,T,HCNE, 1 , 0, L AB A , I T II LE, 1 . C , 5 . 00 , 5 , 1 , 2 , 2 , 3 , 1 0 , 0 , L A ST ) 
CALL DRAW(N, T , HTUO, 2 , 0 , LABT! , ITITLE, 1.0,5.00,5, 1,2, 2, 8, 10,0, LA ST) 
CALL DRAW(N,T,HTHRt,3,0,LAPC, IT I TLE , 1 . 0 , 5 . 0 , 5 , 1 , 2 , 2 , 8 , 1 C , 0 , LA ST ) 
C THIS SECTION COMPUTES SIMPLE Di FF 
READ TAPE 2, HONE 
READ TAPE 2, HTKO . 

READ TAPE 2, HTHRE 
REWIND 2 
DO 15 L=l,900 
H(L)=HONE{L) 

HONE(L) =H(L)-HTWC(L 1 
STORE=HTWO( L ) 

HTKO(L)=H(L)-HrHRE(L) 

15 HTHRE(L)=STORE-HTHRE(L) 

DO 16 L=1 , 12 

16 , ITITLE! L)=8h 

ITITLE! 1 >=8HANDERSQN 
ITITLE!2)=8H BOX 263 
ITITLE!7)=8HFIRST DI 
lTITLE!8)=eHFFERENCE 
ITITLE!9)=8H IN MAG 
ITITLE! 10) = 8HFIELD 
ITITLE! 1 1 )=SHT IN HRS 
ITITLE! 12)=8H H GAMMA 
LABA=4HDEL1 



LABB=4HDEL2 

LABC=4HDEL3 

CALL DRAW!N,T,HONE, 1 , 0, LAB A , ITITLE, 1.0,5. 00, 5, 1,2, 2, 8, 10,0, LAST) 
CALL DRAW !N, T,HTWO, 2,0, LABP , ITITLE, 1.0,5.00,5, 1,2, 2, 8, 10,0, LAST) 
CALL DRAW ! N, T, HTHRE, 3, 0,L ABC, IT I TLE , 1 . 0 , 5 . C , 5, 1 , 2 , 2 , 8 , 1 0 , 0 , L AST ) 
AUTOCORRELATION FUNCTION 75LAGS 
NLAGS=75 
MPTS= N-NLAGS 
PTS= MPTS 



9 



1 1 



XL = 0. 

DO 9 J= 1 , 100 
B I N 1 ! J ) =0 . 
eiN2! J)=0. 

8IN3! J)=0. 

X! J)=XL 
XL = XL+1. 

READ TAPE 2, HONE 
READ TAPE 2,HTW0 
READ TAPE 2, HTHRE 
REWIND 2 
DO 11 J=1,N 

HONEIJ) = HONEIJ)- HONEAV 
HTWOvJ) = HTWC!J) - HTWCAV 
HTHRE! J)=HTHRE! J)-HTHREAV 
DO 12 J=1,NLAGS 
DO 14 K=1,MPTS 
L=K-1 



r> T •> > *« ' « * 



r.tA V «-l lO •> l.CL ^ 1 » 






O T t ^ l_% 



bli\JtJ)= blN^(J) + HTWO(K)<HTWO(L+J ) 

14 QIi\3( J)=RIN?( J) + HTHR£{K)^HTHRL(L+J ) 

BIN1 ( J)=BIK1 ( J) /PTS 
BI,\2(J)= SIN2(J)/PTS 
12 BIN5(J)= BIN3(J)/PTS 
C SAVE AUTOCGRRELATION FUNCTION 
DO 105 J=l,100 

H(J)=eiN1(J) APPENDIX III 

H( J+100)=BIN2( J) 

105 H( J+200)=BIN3( J) 

00 100 J=l, 12 

100 ITITL£(J)=FH 
ITITLEI 1 )=ChANOERSCN 
ITITLE(2)=oH BOX 263 
ITITLE(7»=8HAUTCC0RR 
IT1TLE(8)=PHELAT10N 
ITITLEI 9) =3HFUNCTI0N 
ITITLEI 10)=0H Y IN PR 
ITITLEI 1 1 )=8H00UCTS X 
ITITLEI 12)=8H IN LAGS 
LABA=4H HI 

LABD=4H H2 
LABC=4H Hi 

CALL DRAW(NLAGS,X,BIN1t 1,0,LABA, I TITLE, 0,0, 0,0, 0,0, 8, 8,0, LAST) 
CALL DRAW(KLAGS,X,BIN2, 2,0,LABB, ITITLE,C,0,0,0,0,0,8,8,C,LAST) 
CALL DRAWINLAGS ,X, BIN3, 3,0,LABC, I TI T L E , 0 , 0 , 0 , 0 , 0 , 0 , 8 , 8 , 0 , L AST ) 
C CROSS CORRELATION FUNCTION 75 LAGS 
DO 101 J=1,100 
BIN1 (J)=0. 

BIN2I J)=0, 

101 PIN3(J)=0. 

DO 102 J=1,NLAGS 
DO 103 K=1,MPTS 
L= K-1 

BIN1 ( J)=PIN1 ( J )+ HTWOIK )*H0N5(L+J ) 

BIN2 ( J) =8IN2 ( J ) +HTWO( K) «HTHRE(L+J) 

103 BIN3( J)=BIN3(J)+ HONE ( K ) «HTHRE ( L + J) 

BIN1 ( J)=BIN1 ( J)/PTS 

BIN2(J)= BIN2(J)/PTS 

102 BIN5(J)= B1N3(J)/PTS 
DO 104 J=l,12 

104 ITITLE(J)=SH 
ITITLEI 1 )=8HANDERS0N 
1TITLE(2) =8H BOX 263 
ITITLE(7) =SHCRGSCORR 
ITITLE(8) =8HELATIGN 
ITITL£(9) =OHFUNCTION 

1 TITLE! 10)=8H Y IN PR 
ITITLE! n )=8H0CUCTS X 
ITITLE(12)=8F IN LAGS 
LABA = 4H H21 

LABB = 4H H23 
LABC = 4H HI 3 

CALL DRAWINLAGS ,X, 8 INI, 1,0, LABA, ITI TL E , 0 , 0 , 0 , 0 , 0 , 0 , 8 , 8 , 0 , L AST ) 
CALL DRAWINLAGS ,X,BIN2, 2,0, LABB, ITI TLE, 0 , 0 , 0 , 0 , 0 , 0 , 8 , 8 , 0 , L AST ) 
CALL DRAWINLAGS,X,BIN3, 3,0, LABC, ITI TL E , 0 , 0 , 0 , 0 , 0 , 0 , 8 , 8 , 0 , L AST ) 
C COHERENCY COMPUTATION 75 LAGS 



107 

108 



109 

no 



1 1 1 
106 



DO 106 J = 1,NLAGS 

BINl IJ) = IBIN1IJ)«*2)/IH{J + 100)«HIJ) ) 

IFIBINIIJ)) 107,108,108 
BIN1 I J)=-BIN1 IJ) 

BINl I J) = SORTFI BIN! IJ ) ) 

BIN2IJ) = I BIN2IJ )«*j 2)/I HIJ+i00)»HIJ+200) ) 

IFIBIN2IJ)) 109,110,110 
BIN2 ( J) =-BIN2 I J ) 

DIN2IJ)= SORTFI BIN2IJ ) ) 

BIN3IJ)=I BIN3IJ ) e«2) / (HI J )^H(J+200) ) 

IFI6IN3IJ)) 111,106,106 
BIN3(J)=-EIN3IJ) 

BIN3(J)= SQRTFI BIN3IJ ) ) 

ITITLE(7)=eHC0H£RENC 
ITI TLEI 8) =8HE 
ITITLEI 10)=8H Y IN CO 
ITITLEI 1 1 )=8hH X IN L 
ITITLEI 12 )=8HAGS 

CALL DR AW (NL AGS, X, BINl, 1,0, LABA, ITITLE,0,.2,0,0,0,0,8,8,0,LAST) 
CALL DRAW(NLAGS,X,8IN2, 2,0, LABB, ITITLE,0,.2,0,0,0,0,8,8,0,LAST) 
CALL DRAWINLAGS ,X,B IN3, 3, 0, LABC, I title, 0, .2, 0,0, 0,0, 8, 8,0, LA ST) 
CONTINUE 

REWIND 3 IH-5 

STOP 



END 

MACHINE READ! 1 ARG,2ARG, 3ARG,4ARG) 



TYPEWRITER MESSAGE CCOE TABLE 



/irroivi-ixA XXX 



C0N(T6 =4504041 22050220AP, 
« TB =121203 ! 2C0000C0Ce, 

* 110=1550 121 4C125042CP, 
CON (T12=450404 1220 30220 4E» 

* T1U=141 1 2OOOOOOOOOOC0, 
140=4504040 1501520041', 

« T42=C4140622141630016, 



T7 =11200613010504208, 
T9 =45040412203022048, 
T 11 =1 21 2031 2C00CC00CB ) 
T 13=20062204032604268, 
T15=4504C45112i40120B, 
T4 1 =5406 140 10406030 18, 
T 43 =20224200000000008) 



ft TAPE UNIT ASSIGNMENT TABLE 

ft ASSIGNMENTS FOR CHANNEL 3/4 LIBRARY TAPE 



CONIKO =0 
ft K3 =3320308 

ft K6 =5520208 



K1 =3320^00, 
K4 =3320408, 
K7 =5520508, 



K2 =3320206, 
K5 =5520108, 
K8 =55^0408) 



ft ASSORTED CONSTANTS 

C0NIR4 =73737373737373758, R 1 6=00700000000000008 ) 



ft TAPE 


READING ROUTINE 






IX 


SLJ(K) 


SLJ( L+5) 


. EXIT/ENTRY 


1 ARG 


ZRU(O) 


ZRO(O) 


. TAPE UNIT NUMBER ARG. 


2ARG 


ZRO(O) 


ZRO(O) 


. INITIAL ADDRESS 


3ARG 


ZRO(O) 


ZRO( 0) 


. TERMINAL ADDRESS 


4ARG 


ZRO(O) 


ZRC( 0) 


. MC05, 1 = BINARY, 2 = BC 




LDA7 ( 1 ARG) 


INA(-nB) 






SIUI ( 19X) 


AJP3( L+3 ) 


I CHECK FOR TU ASSIGNMENT 




ENA ( T40 ) 


SLJ4( G) 


. GO PRINT ERROR 




SLJ( P) 


ZRO( 0) 






INA{ IIB) 


STA(KO ) 






LILl (KO)- 


LDA1 (KO ) 






ADD7I4ARG) 


SAU( 4X ) 


I SET TU CONTROL CODES 




SCL(777B) 


SAL( 4X ) , 


. INITIALIZE 




SAL ( 5X) 


SAL( 1 (X) 






INA(3) 


SAU( BX ) 






INA(2) 


SAU( 9X ) 






INA( 1 ) 


SAu ( n X ) 






INA< 1 ) 


SAU (7X ) 






LRS (IS) 


ENA( 0) 






LLS(3) 


SAL( 2X ) 






SAL(3X) 


SAU ( 12X ) 






LLS(39) 


LDQ( R1 Z ) 






STA( KO) 


LDA( 5X ) 






SSU( KO) 


STA( 5X ) 






LDA7(2ARG) 


SAU( 5X ) 






LDA(6X) 


SSU(K0 ) 






STA(6X) 


ENI ( 0) 




2X 


LDA7(3ARG) 


SAL( K) 






SAU( 6X) 


SIL2( 19X) 




3X 


ENI2(4) 


L I L 1 ( N ) 




4X 


EXF (N) 


EXF7(N ) 


I SELECT UNIT 


5X 


EXF ( N) 


EXF7(N ) 


. ACTIVATE CHANNEL 


6X 


EXF(N) 


LDA7(2ARG) 






SUh ( R4) 


LLS(40 > 






INI 1 (-1 ) 


LDAl ( Z ) 






INI 1 ( 1 ) 


SUB(RA ) 




7X 


EXF7IN) 


SLJ ( 13X) 


.* SENSE ENDFILE 


8X 


EXF7(N) 


SLJ( L + 4 ) 


. SENSE PARITY ERROR 


9X 


EXF7(N) 


SLJ( lOX) 


. SENSE LENGTH ERROR 




0JP1 ( 19X) 


AJP( 3X ) 


. SENSE FOR PAD RECORD 




SLJ{ 19X) 


ZRC(O) 


. EXIT 


J 


IJP2( L+2) 


SLJ( 14X) 


. SENSE FIVE PARITY ERRORS 


lOX 


IJP2( L+ 1 ) 


SLJ ( 1 3X ) 


, SENSE FIVE LENGTH ERRORS 




AJP (3X) 


CJP( 3X ) 


. SENSE FOR BAD RECORD 


1 1 X 


EXF (N) 


EXF7(N) 


. BACKSPACE 


12X 


L I L T ( N ) 


SLJ ( 4X ) 


. RETURN TO REREAD 


13X 


£NA(T6 ) 


SLJ ( 17X > 


- SET TO INDICATE TROUBLE 


14X 


ENA( T9) 


SLJ ( 17X ) 




15X 


ENA ( T12 ) 


SL J( IBX) 




17X 


SLJ4( G) 


ZRO(O) III-6 . TROUBLE EXIT 




AJP ( 1 9X ) 


AJP3(5X) 


- . SENSE ACTION TO TAKE 




SLJ ( nx ) 


ZRO(O) 




18X 


SLJ4.; g; 


ZRO( 0) 


I TROUBLE EXIT 




AJP ( 3X) 


AJP3(3X) 


. SENSE ACTION TO TAKE 




SLJ(llX) 


ZRO(O) 




1 9X 


ENI 1 (N) 


EN I 2 ( N ) 


.* retStore indexes 




SLJ (IX) 


ZRO( 0) 





PROGRAM PROMAG APPENDI;C III 

DIMENSION NDATA(864) ,DATE( 100) ,H( 1 7 2 8 ) , HON E ( 900 ) , I TITLE (12), 
1HTW0(900) ,HTHRE(900) ,T( 900) 

DO 5 J = 1,14 
READ 5,DATE(J) 

5 F0RMAT(016) 

CONIMASK = 0000000077777777B) 

ENA(NDATA) , INA( 1 ) . 

STA(INIT).INA(864). 

STAI ITERM) ,£NI (0) . 

6 CALL READ(.^, INIT, ITERM, 1 ) 

00 70 I = 1,11,2 

ADATE = DATE ( I ) 

BOATE = DATE (1 + 1 ) 

LDA( ADATE) ,SUB(NDATA +1 ),AJP1 (6), 

LDA(BDATE) ,SUB(NDATA +2),AJP1 (6). 

NBLOC = 0 

30 CALL READ(3, INIT, ITERM, 1 ) 

ND = NDATA(3) 

LDA( ND) ,AJPO ( 10 ) . 

K=0 

DO 40 J=1,864 
K=K+1 

INDATA = NDATA(J) 

LDA( INDATA) , LDQ ( MASK ) , ARS ( 24 ) , STL ( I NLI ST ) . 

H(K)=INLIST 

K=K+1 

LDA( INDATA) , LDQ (MASK) , S TL ( I NL 1ST ) . 

SLSl (40) . 

40 H(K)=INLIST 
AVE=0. 

DO 50 J =1,1 728 

H( J)=( 1024./H( J) ) *2348400.0 

50 AVE=H(J)+AVE 
AVE=AVE/1728. 

PRINT 26, AVE 

C ZAPP FILTER 

NERR = 0 

AVE = 0. 

HBAR = 50950.0 
DO 51 J=l, 1728 
DEV = ABSF(H( J)-HBAR) 

IF(0EV-200.) 51,52,52 

52 NERR = NERR + 1 
IF(J-3) 53,53,54 

53 H(J) = H(J+3) 

GO TO 51 

54 H(J) = H(J-3) 

51 AVE = AVE + H(J) 

AVE = AVE/1728. 

PRINT 26, AVE 
PRINT 27, NERR 
NBLOC = NBLOC +1 
WRITE TAPE 2, H 
GO TO 30 

10 REWIND 2 

NPTS = ( 1728*NBLOC)/2700 
PRINT 27, NPTS 
READ TAPE 2, H 
K = 1 

TIME = -6.5/900. 

M=1 

C THIS SECTION COMPUTES H VS TIME 
DO 60 L=l,900 

H0NE(L)=(H(K)+H(K+3)+H( K+6) )/3. 

HTWC(L )= (H(K+1 )+H(K+4)+H(K+7) )/3. 

LCOUNT=L 

HTHRE( L )=(H{ K+2 )+H(K+5 ) +H( K+8 ) ) /3. 

TIME = TIME + 6.5/900. 

T(L) = TIME 
K = K + 3*NPTS 
MAX = K+8 

IF(MAX-- 1728) 60,60,80 
80 K=1 

IF(M-NBLOC) 85,95,95 ttt 7 

85 READ TAPE 2,H 
M=M+1 

60 CONTINUE 
95 PRINT 97,LC0UNT 
PRINT 97 -M 

i*t 16 / 100 . 

REWIND 2 

- WRITE TAPE 2, HONE 
WRITE TAPE 2, HTWO 

i.ioTTc TArirr o tjTim c . - . 



r luis t: M V - u . 

HTW0AV=0. 

HTHREAV=0. 

DO U8 L=1,900 

HONEAV=HONEAV + HONE ( L ) APPENDIX III 

HTWCAV=HTWOAV+HTWO(L ) 

48 HTHREAV = HTHREAV+HTHRE (L ) 

H0NEAV=H0NEAV/900. 

HTW0AV=HTW0AV/900. 

HTHREAV=HTHREAV/900. 

PRINT 49tH0NEAVtHTW0AV, HTHREAV 

49 FORMAT! 3E20. 7) 

PRINT 27.NBL0C 

27 FORMAT! 17) 

HONEAV= IHONEAV +HTWOAV +HTHREAV)/3. 

PRINT 26,HONEAV 
26 FORMAT ! E20. 7 ) 

00 90 L=1 ,900 

HONE!L) = HONEIL) - HONEAV 
HTWC(L) = HTWO!L) - HONEAV 
90 HTHREIL) = HTHREIL) - HCNEAV 
DO 1 L=l,12 
1 ITITLE(L) = 8H 

ITITLE! 1 )=8HAN0ERS0N 
ITITLE!2)=8H BOX 263 
ITITLE!7)=8H EARTHS 
ITITLE!8)=8HMAGNETIC 
IT1TL£!9)=8H FIELD V 
ITITLE! 10)=8HS TIME 
ITITLE! 1 I ) = 8HT IN HRS 
ITITLE ! 12 )=8H H GAMMA 
LA8A=4H HI 
LABB=4H H2 
LABC=4H H3 

CALL DR AW! 900, T, HONE, 1, 0,LABA, ITITLE, 1.0, 10. 0,5, 1,2, 2, 8, 10,0, LA ST) 
CALL DRAW1900, T,HTW0,2, O.LABB, ITITLE, 1.0, 10. 0,5, 1,2, 2, 8, 10,0, LA ST) 
CALL DRAW(900,T,HTHRE,3,0,LABC, ITITLE, 1.0, 10. ,5, 1,2, 2, 8, 10,0, LA ST) 
CALL SATGRAF (0, 900,T,H0NE,0,0,0,0) 

CALL SATGRAF (0,900, T,HT WO, 0,0, 0,0) 

CALL SATGRAF !0, 900, T,HTHRE, 0,0, 0,0) 

C THIS SECTION COMPUTES SIMPLE DIFF 
READ TAPE 2, HONE 
READ TAPE 2, HTWO 
READ TAPE 2, HTHRE 
REWIND 2 
DO 1 5 L=1 ,900 
H(L)=HONE(L) 

HONE!L)=H(L)-HTWO(L) 

STCRE=HTWO(L ) 

HTWO!L)=H(L)-HTHRE(L) 

15 HTHRE(L)=STORE-HTHRE(L) 

DO 16 L=l,12 

16 IT1TLE!L)=8H 
ITITLE! 1 )=8HANDERS0N 
ITITLE!2)=8H BOX 263 
ITITLE!7)=8HFIRST DI 
ITITLE!8)=8HFFERENCE 
ITITLE(9)=8H IN MAG 
ITITLE! 10)=8HFIELD 
ITITLE! 1 1 )=8HT IN HRS 
ITITLE! 12)=8H H GAMMA 
LABA=4HDEL1 
LABB=4HDEL2 
LABC=4HDEL3 

CALL DRAW! 900, T , HONE, 1,0,LABA, ITITLE, 1.0, 2. 00, 5, 1,2, 2, 8, 10,0, LA ST) 
CALL DRAW ! 900, T, HTWO, 2, C,LABB, ITITLE, 1.0, 2. 00, 5, 1,2, 2, 8, 10,0, LA ST) 
CALL DRAW!900, T, HTHRE, 3,0, L ABC, ITITLE, 1.0, 2. 0,5, 1,2, 2, 8, 10,0, LAST) 
CALL SATGRAF !0 , 900 , T , HONE , 0 , 0 , 0 , 0 ) 

CALL SATGRAF ! 0, 900, T, HTWO, 0,0, 0,0) 

CALL S ATGRAF ! 0, 900, T, HTHRE, 0,0, 0,0) 

70 CONTINUE 
REWIND 3 
STOP 
END 

MACHINE READ! 1 ARG,2ARG, 3ARG,4ARG) 



* NOTE G = ADDRESS OF ENTRY TO TYPEWRITER SENSE ROUTINE 

* P = ADDRESS OF PROGRAM CONTROL ROUTINE JUMP ENTRY 

* THIS ROUTINE IS FOR PROGRAM CALL 

L0C!G*31, P = 10) IJl-S 

* TYPEWRITER MESSAGE CODE TABLE 





m 








• w 



I fc- I -r . 



CON (T12=45040U1 22030220 4B, 

* T 14=141 1 200000000000B, 

* T40=4504040130152004B, 

» 142 = 041406221416300 IB, 






113=20062204032604266, 
115=45040431 12 1401 20B, 
T41=3406140104060301B, 
T 4 3=20224200000000006) 



» TAPE UNIT ASSIGNMENT TABLE 

* ASSIGNMENTS FOR CHANNEL 3/4 LIBRARY TAPE 



CONIKO 
» K3 

► K6 



=0 , 
=3320306, 
=5520206, 



K1 

K4 

K7 



=3320106, 

=3320408, 

=5520306, 



K2 

K5 

K8 



=3320206, 
=5520106, 
=5520406 ) 



* ASSORTED CONSTANTS 

CONIRU =73737373737373736, R16=0070000000000000B ) 

* TAPE READING ROUTINE 



IX 


SLJ (N) 


SLJ (L+5) 


. EXIT/ENTRY 


lARG 


ZRO(O) 


ZRO(O) 


. TAPE UNIT NUMBER ARG. 


2ARG 


ZRO(O) 


ZROIO) 


. INITIAL ADDRESS 


3ARG 


ZRO(O) 


ZRO(O) 


, TERMINAL ADDRESS 


4ARG 


ZRO(O) 


ZRO(O) 


. MODE, 1 = BINARY, 2 = DC 




LDA7{1 ARG) 


INA(-1 IB) 


• 




SIUl ( 19X) 


AJP3IL+3) 


. CHECK FOR TU ASSIGNMENT 




ENA(T40) 
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